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Abstract. We consider the Hartree-Fock approximation of Quantum Electro- 
dynamics, with the exchange term neglected. We prove that the probability 
of static electron-positron pair creation for the Dirac vacuum polarized by an 
external field of strength Z behaves as 1 — exp(— kZ'^/^) for Z large enough. 
Our method involves two steps. First we estimate the vacuum expectation of 
general quasi-free states in terms of their total number of particles, which can 
be of general interest. Then we study the asymptotics of the Hartree-Fock 
energy when Z — >■ -\-oa which gives the expected bounds. 
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Introduction 

In 1930, Dirac jB] suggested the idea of identifying the vacuum with a sea of 
virtual electrons with negative kinetic energy. His theory implies that when a suf- 
ficiently strong source of energy is provided to the vacuum, some virtual electrons 
are excited into real electrons, leaving "holes" in the Dirac sea. These holes can be 
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interpreted as positrons, the anti-particles of the electrons, which were experimen- 
tally observed in 1933 by Anderson [1]. The extraction of an electron from the Dirac 
sea is usually called electron-positron pair creation. Sauter [26j . and Heisenberg- 
Euler ^ISj considered the possibility that an external electromagnetic field could 
excite the Dirac sea to create those pairs. Schwinger [28] then computed the prob- 
ability of dynamical pair creation by a constant, uniform, external electric field in 
the framework of Quantum Electrodynamics (QED). The specific phenomenon of 
pair production triggered by an external, non-quantized field is thus labeled the 
Schwinger effect. It is remarkable that this effect is different from the absorption 
of photons by the vacuum, which is another possible source for pair creation. The 
Schwinger effect is based on the fact that the vacuum acts as a polarizable medium 
which can decay into electron-positron pairs when excited by a sufficiently strong 
electric field. Although the modern formulation of QED no longer describes the 
vacuum as a sea of virtual particles, Dirac's theory is still valid in the mean-field 
approximation [16J. 




Experimentally, pair creation in electric fields has not been observed yet because 
it is only non-negligible in a very strong field. However, recent progress in laser 
physics have permitted to create very strong fields, making the observation of the 
Schwinger effect possible in the near future [91 [301 [S] . 

One has to distinguish between dynamical and static pair creation. Dynamical 
pair creation consists in studying the time evolution of the vacuum state when an 
external field is progressively turned on, so that a pair consisting of a scattering 
electron and a corresponding hole in the Dirac sea is created. The external field is 
then progressively switched off, and one has to check if the pair still exists when the 
field is completely turned off. Static pair creation, on the other hand, consists in 
the study of the absolute ground state (the polarized vacuum) of the Hamiltonian 
in an external field. Therefore, it is a time-independent process. In this context, 
the vacuum with an additional particle is energetically more favorable than the 
vacuum without particle. Static pair creation is easier to study than dynamical 
pair creation, but it is also a bit less relevant from the physical point of view. 

When the interactions between particles are neglected (the so-called linear case) , 
static pair creation was mathematically studied by Klaus and Scharf |20| . They 
proved that the probability of pair creation becomes 1 when the strength of the 
positive external field sufficiently increases such that an eigenvalue of the Hamil- 
tonian of the system crosses zero. In the linear case, dynamical pair creation is a 
very involved phenomenon, whose properties were mathematically understood very 
recently. Nenciu [22, 23 proved that there is a discontinuity in the probability to 
create pairs as the strength of a specific external field increases, in the adiabatic 
limit. Later on, Pickl and Diirr JH [24] proved that the probability of pair creation 
tends to 1 in the adiabatic limit, for general over-critical external fields, by care- 
fully studying the resonances created by the eigenvalues diving into the essential 
spectrum of the Hamiltonian of the system. 

This article is devoted to the mathematical study of static pair creation in a non- 
linear model describing the polarized vacuum, taking into account the interactions 
between particles. This model was first proposed by Chaix and Iracane [6] in 1989 
and it has recently been given a solid mathematical ground in a series of papers by 
Gravejat, Hainzl, Lewin, Sere, and Solovej [13l [Mj [171 [151 [10] . As in those papers, 
the main difficulty of our work is the nonlinearity of the model. The more involved 
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study of dynamical pair creation for the same model will be the subject of future 
work. 

In the considered model, the polarized vacuum in a potential generated by a 
density of charge Zv is described by an operator Pz on _L-^(M^,C*) (a density 
matrix). This operator is a solution to the nonlinear equation 

\ = Xr-oo.o](^z) + '^ 

\Dz = i^" + a (pp,_ 1 - Zi.) * I • ' 

where I?" :— —ia. • V + /? is the (free) Dirac operator and (5 is any self-adjoint 
operator such that < (5 ^ 1 with rank(5) C ker(£>2). Discarding the operator 5, 
we see that Pz is the ground state in the grand canonical ensemble of a dressed 
Dirac operator with density of charge Zv perturbed by the density pp^^i of the 
vacuum. The polarized vacuum therefore interacts with itself. 

We shall consider the operator Pz in the limit Z — > +oo. Of our particular inter- 
est is the probability that pairs are generated, which is a nonlinear function of Pz 
(see Section 11.11 below) . The usual picture [20l |27l [12] is that if the first eigenvalue 
Xi{Z) oi Dz is negative, as showed in Figure [U then the vacuum becomes charged 
and the probability of creating at least one pair is 1. In the linear case, Hainzl [12] 
showed that the charge of the vacuum in the external density Zv is exactly the 
number of eigenvalues (counted with multiplicity) of the operator — tZv-k | • 
crossing when we increase t from to 1. However, because of the nonlinearity of 
the model we study, detecting for which values of Z the first eigenvalue will cross 

is very difficult. However, the probability of pair creation can be very close to 1 
without any crossing, as we will explain in Section 11.11 

Ai(Z) 

-1 1 

Figure 1. Spectrum of the mean- field one body Hamiltonian Dz- 

More precisely, we prove that the probability of static pair creation behaves as 

1 — exp(— kZ^/'^) (see Theorem [T]), where Z is the charge of a nucleus put in the 
vacuum, and k is a constant depending on different parameters of the model such as 
the cut-off or the shape of the nucleus. The proof relies on the large- Z asymptotics 
of the polarized vacuum energy, which is obtained by using an appropriate trial 
state. This implies that the average number of particles of the polarized vacuum is 
of order Z^/^ . We then use general estimates showing that the probability to create 
pairs for a quasi-free quantum state is bigger than 1 — exp(— k'A^), where N is the 
average number of particles of the quantum state and n' is a universal constant. 
Since for the polarized vacuum TV ~ Z"^/^ , the result follows. 

The paper is organized as follows. In Section 1 we introduce the Bogoliubov- 
Dirac-Fock model and we state our main result. In Section 2, we prove the general 
estimates on quasi-free states on Fock space, which are of independent interest. In 
the end of Section 2 we come back to our particular setting. In Section 3, we study 
the large-Z asymptotics of the polarized vacuum energy. Finally, in Section 4 we 
prove Theorem [T] using the tools developed in Section 2 and 3. In Appendix A, we 
recall some properties of product states, which are used in Section 2. 



4 



JULIEN SABIN 



Acknowledgments. I sincerely thank Mathieu Lewin for his precious guidance 
and constant help. I also acknowledge support from the ERC MNIQS-258023 and 
from the ANR "NoNAP" (ANR-IO-BLAN 0101) of the French ministry of research. 

1. Estimate on the probability to create pairs 

1.1. Probability to create a pair. The first quantity to define is the probability 
to create a pair. Let S)+^S)- be (separable) Hilbert spaces, representing the one 
particle (resp. anti-particle) space. The natural space to describe a system with an 
arbitrary number of particles/anti-particles is the Fock space 

with the usual notation F{S)) := ©at^o f^i -ft for any Hilbert space and with 
the convention A^i^ := C We also define the vacuum state := r2+ ® r2_ G Jvi 
where Vt± := 1 ® ffi • • • € J^{S)±). Recall that a state over can be definedlj 
as a positive linear functional lo : B{Tq) — > C with a;(Idjro) — 1, where B{To) is 
the set of all bounded linear operators on J^q. Notice that any normalized ^ € J^q 
defines a state (called pure state) by the formula u)^{A) = {i/j , AiJ;) jr^ , where 
{■,-}j^o is the usual inner product on J"o- Following [24l Corollary 4.1], [311 Eq. 
(10.154)], and [131 Section 2], we define the probability p{uj) for a state w to create 
a particle/ anti-particle pair by 

(1.1) \p{uj) ■.= l-u;m{n\), 

where G B{J-o) is the orthogonal projection on Cil. For a pure state tp = 

V'o,offiV'o,i®V'i,offi • • • G J-Q, we have p(a;^) = 1 — |^o,o| .Therefore, _p(w^) = if and 
only if ■0 = (the vacuum has probability zero to create pairs), while p{uj^) = 1 
if and only if ^/;o,o = 0- In the latter case, notice that ■0 does not litterally contain 
pairs, in the sense that its number of particles may not be equal to its number of 
anti-particles. This definition merely measures the probability that a state contains 
real particles/anti-particles. 

Typically, f2 represents the free (or bare) vacuum and we want to measure the 
probability of a perturbation i7' of f2, representing the polarized (or dressed) vac- 
uum in the presence of an external electric field, to have pairs. Assuming that 
is a pure quasi-free state, we have the well-known formula (see e.g. [3T1 Theorem 
10.6], [51, Theorem 2.2], or Theorem 5]) 

TV M 

(1.2) n' = l[ ^— f n n ^0i9rn) [] (1 + >^^<iv^K{u^)) O, 

i^l V 1 + n=l m=l i^l 

where Oq (resp. 5q) is the free particle (resp. anti-particle) creation operator, 
{fn)n U {vi)i (resp. {gm)m U {ui)i) is an orthonormal set for Jo+ (resp. io-), and 
(A,), e e{R+). 

From the formula (|1.2|) . we see that p(ajsi') = 1 as soon as > or M > 0. 
Moreover, in this case real particles in the states (/»)« and real anti-particles in the 
states {g,n)m have been created. In the linear case, Klaus and Scharf [20] proved 
that N, M 7^ if the external field is strong enough. However, there can be a high 



For convenience, we will later define a state as a linear form on the tensor product of CAR 
algebras Cgi — ] C B(J-'o), which in our context does not change anything since all the 

states we consider are normal. 
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probability to create pairs even when N = AI — 0. Indeed, since in this case we 
have 

one sees that ^(wq/) is close to 1 when the are large enough. One simple condition 
is that J2i A? (1 + is large enough, by the inequality 



- + A 



cxp 



+ A 



Note that this is indeed (half) the average total number of particle of the state fl' 
(number of particle + number of anti-particle), 



A? 



where Af is the usual number operator on J-q (see formula (j2.15p ). Hence p(ll)(i>) 
is close to 1 when wq' {Af) is large enough. While the non- vanishing N, M case 
can be interpreted as the creation of real particles, this second explanation for an 
increasing p(wsi') can be interpreted as a "virtual pair creation". In this article, we 
study an analog of "virtual pair creation" for more general states than those given 
by formula pT2)) . 



1.2. Static pair creation in the reduced BDF approximation. For nonin- 
teracting electrons in an external field V , the polarized vacuum fi' is the unique 
Hartree-Fock state whose density matrix is [20l [12] 

^-X(-oo.0](^°+^). 

In this article, we will rather use the reduced Bogoliubov-Dirac-Fock approximation, 
a non-linear model enabling to describe an interacting vacuum in which y is a 
function of P itself. It was introduced by Hainzl, Lewin, Sere and Solovej in a 
series of articles [131 IHl Ull HI] after the pioneering work of Chaix, Iracane, and 
Lions [H [7] . We will now briefly recall the model and the results needed for our 
study. 

In units where m — c — h — \, the reduced Bogoliubov-Dirac-Fock (rBDF) 
energy functional is the (formal) difference between the energy of the state P and 
that of the free vacuum P^ = X(-oo,o] (-C"), with the exchange term dropped. It 
depends only on the variable Q = P — P^l , 

(1.3) ■■^TMD''Q)~aD{pQ,:.) + ^D{pQ,pQ). 

Here, a > is the coupling constant and j/ : R'^ — M is the external charge density 
belonging to the Coulomb space 

C:=lfe.y'(R'): I %ll!dfc<+oo|. 



endowed with the inner product D{pi, P2) = J \k\ pi{k)p2{k) (the hat denotes the 
Fourier transforrrQ) . We also use the notation = D{p,p) /"^ for any p <eC. In 



^Our convention is f(k) = (27r)-»/2 f{x)e-^^-'^ dx, Vfc e 1 
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order to define the domain of the rBDF energy functional, let us fix a cut-off A > 
and define the one-particle Hilbert space 

JOA := {/ei'(K^C4), supp7cS(0,A)}. 

The operator 1)° = —icx - V + j3 \s the usual Dirac operator on L^(R'^, C"*), where 
ai, q;2, aa, /3 are the Dirac matrices acting on C^, 

and (cri)i=i.2,3 are the Pauli matrices, 

Ol\ f Q -i\ /lO 



The operator stabilizes i^A, and its restriction to Joa defines a bounded 
operator on j^A, which we still denote by I?". For convenience we introduce :— 
1 — P^. We denote by &p{S)) the Schatten class of all bounded operators A on the 
Hilbert space such that Tr(|A|P) < oo. For any operator Q on i^A and for any 
£,£' e {+, -}, we let Q^^, P°QP°, and we define 

e.^po (i3A) {Q e ©2(iOA), Q++, Q— e ©i(i3A)} ■ 

It is a Banach space endowed with the norm 

IIQIIi.p^ IIQ++llei + HQ— llei + IIQ+-Ile2 + IIQ-+lle2- 

For any Q G ©i.po {^a), we define its generalized trace by 

Tro(g) :=Tr(Q++ + g__), 

and its density pq by Pq{x) := Ttc4.{Q{x,x)) for all x G M.^, where Q{x,y) denotes 
the 4x4 matrix kernel of Q. This density pq is well defined since supp(3(-, •) C 
5(0, A) X B{0,A) implies that Q{-, •) is smooth. Furthermore, it is proved in [151 
Lemma 1] that pq e L^(]R^) n C for any Q G &i po {^a)- We conclude that the 
rBDF energy functional is well-defined on the convex set 

(1.4) /C := {q e 6i,po(ioA), Q^Q*, -P°^Q^1-P^]. 
Notice that the kinetic part of the rBDF energy is well defined since 

Tro(i?°Q)-Tr(|i?0|(Q++-Q__)). 
The variational set JC is the convex huU of {P - P°, P = ^ p* p _ pQ ^ 
&2{^a)}, where P is the density matrix of a pure Hartree-Fock state, which is a 
Hilbert-Schmidt perturbation of the free vacuum P^. The rigorous derivation of 
the rBDF energy functional and the motivation for this functional setting can be 
found in [HI [12. 



For any Z > and v G C, the rBDF energy functional S^^j^p admits global min- 
imizers Qz on JC. Minimizers are not necessarily unique, but they always share the 
same density pz '■= Pqz- Any minimizer Qz satisfies the self-consistent equation 



(1.5) 



Qz - X{-oc.o]iDz) - P"^ + S 
Dz = D"^ + a{pz-Zv)^\■\- 



v^\leYe i5 is a self-adjoint operator such that < (5 ^ 1 and rank(5) C ker(P>2). 
Hence, uniqueness holds if and only if \ev{Dz) — {0}. Notice that since the density 
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pz is unique, the operator Dz is itself unique. Any minimizer of E^i^p on K, is 
interpreted as a generalized one-particle density matrix of a BDF state oj^^^ (see 
Section [2. 6p representing the polarized vacuum in the potential Zv-k \ ■ \^^. When 
there is a unique minimizer Qz, it is a difference of two projectors and hence it 
is the generalized one-particle density matrix of a pure state. We emphasize that 
there is no charge constraint in this minimization problem: In fact, the polarized 
vacuum could (and should) be charged when Z is very large. In this case, one may 
think that an electron-positron pair is created, with the positron sent to infinity. 

We want to estimate p{uj^^^) in terms of Z, and confirm the picture that the 
stronger the field, the more pairs are created. As a consequence, we will fix a non- 
zero density v (interpreted as the shape of the external charge density) and study 
Pz ■— p[ijJ^^c) for Z > large. Our main result is the following. 

Theorem 1. Let a > and A > 0. Let £ C such that J^g^sil + |a;|)|i^(a:)| dx < oo 
and q := J^^ v ^ 0. Then, there exists a constant Zi > and a constant k > 
such that for all Z > Z\ we have 



(1.6) 



. '72/3 



The constant Z\ is equal to Z^j J^g \v\ where Z^ is defined in Equation (j3.8p and 
Z\ depends only on A, a, \q\, and J^j The constant k equals 0.0941248 

C(Jju3 |j/|)^/'^, where C is defined in Equation (|3.9|) below. 



X 



Remark 1.1. The assumption \x\iy £ L^{M.^) allows us to have an explicit estimate. 
If we remove this assumption, we can still prove the weaker result that — 1 as 
Z oo. In the sequel, by rescaling Z if needed, we will assume that 



kl = l. 



If J v = 0, we expect an asymptotics lower than Z^/"^, but we are unable to prove 
it. 

Remark 1.2. We will see that Zi ^ const, x a^"^/^ as a — > 0. 

Theorem [1] says that in a very strong field, Z ^ 1, the probability to create 
at least one electron-positron pair is very close to 1. It is reasonable to think 
that for some sufficiently large Z, the first eigenvalue of Dz crosses in which case 
Pz — 1. However, determining the behaviour of the eigenvalue of Dz as Z increases 
is difficult because of the nonlinearity of the model and because we are in a regime 
far from being perturbative. For all these reasons, the estimate (|1.6p on pz is the 
best we have so far. For very large Z, one expects that many electron-positron pairs 
will be generated. We conjecture that we have indeed uj^^^{Vk) as Z ^ +oo 
for all k G N, where Vk is the orthogonal projector on the fc-particle space in Fock 
space (see Section [5]). This would mean that for large Z, the probability to create at 
least fc- pairs is very close to 1. Our method of proof only gives this result for fc = 0. 
However, if we assume that w^"^ is a pure quasi-free state for all Z large enough 
(which is the case if ^ (j{Dz)), then the conjecture follows from Proposition 12.61 
below. 
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1.3. Strategy of the proof. The proof is separated into two parts. The first one 
consists in estimating the energy of the polarized vacuum, 

(1.7) EiZ) := £:fBDF(Qvac) = inf {S^MQl Q & IC} , 

from above by —cZ^^^. We will also give a lower bound E{Z) > —Z^^^ to show 
that the power 5/3 is optimal, although we only need the upper bound for the 
proof of Theorem [T] From this estimate, we then infer that the average number 
of particles (counted relatively to that of the free vacuum, see Section 12. 6|) in the 
polarized vacuum satisfies 

Tr((Qi''J++-(04''J__)>Z2/^ 

The precise statements of these results and their proofs can be found in Section |31 
In a second part, we prove an estimate on the vacuum expectation w(|r2)(r2|) 
for a quasi- free state w, in terms of its average number of particle a;(A/'). These 
estimates are of independent interest and therefore we also provide several other 
estimates for the distribution of quasi-free states in the /c-particle spaces. These 
results are contained in Section [2l Finally, in Section |4l we combine the two parts 
and prove Theorem [1] 

2. On the distribution of quasi-free states in the fc-PARTICLE SPACES 

In this section, we consider general quasi-free states. Only in Section 12.61 we 
come back to our particular situation of pair creation. We start by introducing the 
notation used throughout this section. 

2.1. Notation. Let (i5, (•,•)) be a complex, separable Hilbert space whose inner 
product (•, •) is linear in the second argument. We also need an anti-linear operator 
J : ^ such that J* J = Idj^, where (.ft, (•, •) j^) is another complex Hilbert spac^l. 
Let F := (Bn^o^^ be the associated Fock space with S^^ := A^^}. We stiU denote 
byf2 = l®0©---eJ" the vacuum vector. For fc € N, we denote by Vk & the 
orthogonal projection on S^'' C J-. We recall from Section [TTT] that B{F) is the space 
of all linear bounded operators on F. Let 2l[jo] be the CAR unital C*-subalgebra 
of B{F) generated by the usual creation (resp. annihilation operators) a* (/) (resp. 
a{f)), ior f £ Sj. We denote by Af the particle number operator on F, 

AA:=0A;ld^. -^a*(/,)a(/,) 

for any orthonormal basis {fi)i^o in Sj. Then Vk — l{AA=fc} for all A: e N. A state 
on 2l[io] is a non-negative linear functional uj : 2t[i3] — > C which is normalized: 
uj{ldjr) = 1. A state w is called normal if there exists a non- negative operator 
G on (sometimes called the density matrix of oj) such that Trjr(G) — 1 and 
uj{A) = Tijr{GA) for all A £ ^[Sj]. Of particular interest are the pure states which 
are normal states with G = \tp) {tp\ ior ip E F with WtpW jr = 1. We define the average 
particle number of uj as 

oj{U) ^c.(a*(/,)a(/,)) e [0,-|-oo]. 



Recall that the adjoint J* of an anti-linear operator is defined as {J*f,g) := {Jg,f)s% for all 
/ G .ft and g £ Sj. Typically, one chooses ^ = S} and J the complex conjugation, or .S = ^i* and 
J(/) = (/, •) 1291 . Here we keep ^ abstract because this will be useful for the construction of BDF 
states in Section IZ6l 
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The one-particle density matrix (1-pdm) 7 of w is the operator defined by 

(g,7/) ■.^u;ia*if)aig)), 

for all f,g E S). It is a self-adjoint operator on S), satisfying ^ 7 ^ 1. In the same 
fashion, we define its pairing matrix a : which is a linear operator on ^ by 

{aJf,g) ■.^uj{a*{f)a*{g)), 

for all /, 5 G io. It satisfies (a J)* = —a J. Moreover, if we define the operator 
r(7, a) on i3 ® J? by block 

(2-^) r(7,«):=(; ), 

then ^ r(7, a) ^ 1, see [U Lemma 2.1]. This last relation implies that 

(2.2) 7^ + aa* 7, 

in the sense of quadratic forms on ^. Notice also that uj{Af) — Tr(7). A state uj is 
called quasi-free if for any operators ei, . . . , e2p which are either a a*(/) or a a{g) 
for any f,g € f), then 10(6162 ■ ■ ■ 62^-1) = for any p ^ 1 and 

(2.3) a;(eie2 . . .e2p) = ^ (-I)^^''''w(e7r(i)e7r(2)) ■ • ■ w(e^(2p_i)e^(2p)), 

7rG(S2p 

where iS2p is the set of permutations of {1, . . . ,2p} which verify 7r(l) < 7r(3) < 
• • • < 7r(2p — 1) and 7r(2j — 1) < 7r(2j) for all 1 ^ j ^ p, and e(7r) is the parity 
of the permutation tt. The relation (|2.3p is called the Wick formula. From this 
definition, we see that a quasi-free state is completely determined by its density 
matrices (7, a). We recall [H Theorem 2.3] 

Proposition 2.1. For any (7, a) such that ^ r(7, a) ^ 1 with additionally 
Tr(7) < -fcx), there exists a unique quasi-free state lj on 2l[i5] with finite number 
of particle such that 7 is its 1-particle density matrix and a its pairing matrix. 
Furthermore, u is normal: there exists G : T ^ T with ^ G ^ 1 and Trjr(G') = 1 
such that uj{A) = Tr^(GA) for all A e 2l[io]. 

Now, we need some terminology, which is not universal in the literature. We call 
Hartree-Fock (HF) states the quasi-free states with Tr(7) < +00 and a ~ 0, because 
when such states are pure, they are usual Slater determinants. Quasi-free states 
with Tr(7) < -\-oo and a 7^ 0, are called Hartree-Fock-Bogoliubov (HFB) states. 
Pure HFB states are particularly simple since they are Bogoliubov rotations of the 
vacuum Q. The aim of this section is to study the distribution of quasi- free states 
in the particle subspaces S)'^, in terms of Tr(7). Our results are different for HF or 
HFB, pure or mixed states. 

2.2. Motivation. Quasi-free states are also called Gaussian states, in particular 
because they can be written as (limits of) Gibbs states of quadratic Hamiltonians 
(i.e. normal states with density matrices e~^^/ Tr(e~^'*), where H is a quadratic 
Hamiltonian). The Gaussian character of quasi- free states is however deeper. In this 
section we will show that the distribution of a quasi-free state w over the different 
S)'', that is (w('Pfe))fc^Oj also has some Gaussian characteristics. More precisely, we 
will provide estimates of the form 

uj(rk)^ck6--''^^^\ 
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This estimate means that a quasi-free state which has a large average number of 
particles uj{J\f) ^ 1 necessarily has an exponentially small vacuum expectation 

oj{\n){n\) ^coe-^'^^^l 

Let us explain the picture in a commutative setting. Let f{x) = Tr^i/^e^'^^"' 
for a e be a Gaussian function such that J^f — 1. Then a — J^xf{x)dx 
is the average position of /, as desbribed in Figure [5] Now if a goes to +oo, 
the whole function moves to infinity and in particular /(O) becomes smaller and 
smaller: /(O) = 7r^^/^e^° . In other words, as the average position of / goes to 
infinity, /(O) goes to zero (and this is true for any f{xo) with xq fixed). We will 
prove a similar fact for quasi-free states. Indeed, for any quasi-free state lo we have 
1 = a;(Id_7r) = J2k^o^('^k), which is the analog of Jg^f{x)dx — 1. We also know 
that a;(A/') = X]fc>o ^^{'^k) is the average number of particle of w; it is the analog 
of dx. We want to prove that when uj{J\f) is large, then the main part 

of CO lives in the high-A; particle spaces, that is a; "follows" its average number of 
particles, as shown in Figure [5] The analog of /(O) in this case is w(|r2)(r2|) and 
we thus want to prove that goes to zero as uj{Af) goes to +oo, for any 

quasi-free state uj. A natural extension of this result would be that uj{Vk„) also 
goes to zero for any fixed k^. 




Figure 2. Analogy between a Gaussian function and a quasi-free state 



We will provide explicit estimates depending on the properties of the quasi-free 
state (pure, mixed, HF or HFB). In the most general case of mixed HFB states, we 
only derive a bound on the vacuum expectation. The following table tells us where 
each case is treated. 





Pure 


Mixed 


HF [a = 0) 


Section O 


HFB [a ^ 0) 


Section [23] 


Section [23] 



In spite of their usefulness, we have not found the following estimates in the 
literature. One main reason is probably that e~^^ does not belong to the CAR 
algebra, hence u}{e'^^^) only makes sense for normal states. 

Remark 2.2. A useful tool for the proofs of the following results is the notion of 
product state. A product state ®iUJi is a state on J-{®iS^i) ~ ®iJ'{S^i) when each 
(jji is a state on F{Sji). While this notion is intuitive, we recall how it is precisely 
defined in Appendix VK\ 
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2.3. Hartree-Fock case. 

Proposition 2.3 (HF case). Let lo he a quasi-free state with Tr(7) < +oo and 
a = 0. Then for any /3 ^ we have 



(2.4) 



(e-^^) = Det^ (1 + (e-f^ 1)7) 



We also have the following estimate 
(2.5) 



for all k ^ ka while uj(Vk) =0 if k < k^, where fco diniker(7 — 1). 

Remark 2.4. This estimate implies that for any fixed fc G N, uj{'Pk) — > as Tr(7) 
+00, which is the expected behaviour. 

Remark 2.5. A more theoretical corollary of (12. 4p is that for any fixed / : N — )■ C 
vanishing at infinity, oj{f{J\f)) goes to zero as Tr(7) goes to +00. This uses the fact 
that the algebra of these fs is generated by the (e~^ )^^o- In the same fashion, one 
can also prove that uj{K) — ^ as Tr(7) — !• +00, for any fixed compact operator K. 

Proof. Since 7 is trace-class it can be diagonalized in an orthonormal basis {fi)i^fi, 
7 = X^i^o For all i, let oji be the unique quasi-free state on J"(C/i) 
having AJdc/i as 1-pdm and Oc/; as its pairing matrix. Then by Proposition IA.2I 
in Appendix 1X1 one has uj — ®iUJi. Moreover, since 

Te-I'^T* = (g) (1 + {e~^ - l)a* (/,)a(/0) , 



where T is the isometry between F{®iS^i) and defined in Appendix lA} we 

have 



(1 + (e-^ - l)a*(/,)a(/.)) - 1] + ^'^'^ 



Detf, (1 + (e 



1)7) 



To prove ([23]), we notice that for aU ,3^0, a;(e"'^-^) = Y^k^n e^^''uj{Vk), and we 
identify the coefiicients of e"'^*'' in HieN + (^"'^ ~ l)Ai). This yields 



uj{Vk) 



icN lei jti 



/cN Vie/ 
#i=k 



^Tr(7) 



ERA. 

/CN 16/ 



-Tr(7) 



E An 



ii,...,ifc6N 

where we used that ^ Ai ^ 1 for all i. 
oj{Vk) = if fc < dimker(7 - 1). 



_ (eTr(7))'- 



Notice from the first equality that 

□ 
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2.4. Pure Hartree-Fock-Bogoliubov case. 

Proposition 2.6 (Pure HFB case). Letuj a quasi-free pure state with Ti'{j) < +oo. 
Then for any /3 ^ we have 



(2.6) 



uie-"-^) = Det^ y^l + (e-2/^ - 1)7. 



We also have the following estimate for all k — kg + 21 with fco := diniker(7 — 1) 
ande>0 



(2.7) 



while uj{Vk) ^Oifk<ko or k ^ ko + 2£ + 1. 

Proof. It is well-known 4, Theorem 2.6] that ui is pure if and only if r(7, a)^ = 
r(7,a), which is equivalent to 7^ + aa* = 7 and [7, a J] = 0. The operator 7 is 
trace-class and aJ is anti-hermitian and Hilbert-Schmidt, hence both 7 and aJ 
are diagonalizable. Since they commute, they are simultaneously diagonalizable. 
Remember that any anti-hermitian can be diagonalized in 1 x 1 blocks corresponding 
to its kernel and 2x2 blocks. Hence there exists a decomposition = (Bi^o^i, 
with dini{S)i) ^ 2 such that 

• For all i, 7 and a J stabilize ^t; 

• If dim(ioi) = 1 then 7|^. — Aildy-,. , aJ|y-,. ~ 0; 

• If dim(jOj) = 2 then 7|j--,^ = o' A ) '"-^l^^. = 
a; e K and a? = Ai — A? 











with 



In particular, uj — ®i^oLJi where uji is the quasi- free state on J-{f)i) with 1-pdm 
7|j5|. =: 7i and pairing matrix aJ^f^.. Let us now prove that for all i 



(2.8) 



,(e-^-^0 = Det,,. ^1 + (e-^P - 1)7, 



where A/^ is the number operator on J^{:P)i). First we consider the case dim{Sji) — 1, 
and let fi G i^i be a normalized vector. Then J-"(i3i) — C® C/i and M = a* {fi)a{fi) 



so that e"''-^' = 1 + (e 



-l)a*(/,)a(/,). Therefore 
a;,(e-^^')-l + (e"^-l)Az 



Since 7^ -I- aa* = 7, we have = if dim(i3i) = 1 hence A^ 
cases we have 



or Ai = 1 . In both 



uj,{e-^^^) = Det,,. ^l + (e-2/^-l)7,. 

Now suppose dim{Sji) = 2 and let {fi,gi) be an orthonormal basis of such that 
in this basis 7|^. and aJ^j-^. have the form given above. Then J-{f)i) = C © C/j © 

Cg^ © C/i A gi and M = a*{fi)a{fi) + a*{gi)a{gi), so that 

e^^^' - (1 + (e-^ - l)a*(/.)a(/,))(l + (e"^ - l)a* (50a(5^)). 
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We deduce that 

co^ie-^^') - l + (e-^-lK(M) + (e-^-l)2^,(a*(/,)a(/,)a*(ft)a(ft)) 

= 1 + (e-^ ~ l)c^^(M) + (e-^ - l)^[u;,ia*{,U)a{f^)Hia*i9iHg^)) 
-Wi(a*(/i)a*(5i))a;j(a(/i)a(5i)) + Wi(a*(/i)a(gi))wi(a(/,)a*(5^))] 

= 1 + 2(e-'3 - 1)A, + [e-^ - ifX, 

= 1 + {e-^^ - 1)X, = Dct^,, .^l + (e-20-l)7„ 

where in the second equaUty we used Wick's relation for Ui. The equahty ()2.6p 
then follows by taking the product of the relations (|2.8p . Putting aside the indices 
i such that Ai = 1 , we obtain 

u;{e-^^)=e-^''"l[{l + {e-'^-l)X.), 

where fco dimker(7 — 1). Identifying the coefficient of e"'^'^ in both sides, as in 
the proof of Proposition [231 we find that ijj{Vk) = for all k < ko, and that 



/CN lel i4I 

/Tr(7) — fcoV Tr(^)^fco e^/"^ /Tr(7)\^ Tr(T) 



for k = ko + 2£ with £ ^ 0. This concludes the proof of Proposition 12.61 □ 

2.5. Mixed Hartree-Fock-Bogoliubov case. In the most general case of a mixed 
HFB state, we cannot apply the same strategy as in the previous cases, i.e. identify 
a; as a product of states living on smaller dimensional spaces. Indeed, 7 and a J can 
have no common stable finite-dimensional subspaces. However, we can still prove 
an estimate on the vacuum expectation. 

Proposition 2.7 (Vacuum Expectation of a Mixed HFB State). Let uj a quasi- free 
state with Tr(7) < +00. Then we have the following estimate 



(2.9) 



where a = max^^o 4e^f>'+7el\ 0.0941248 . . .. 

Remark 2.8. We believe that for any sequence of quasi-free state {ujn)n with finite 
number of particles, we have ujn{e~^^) — >■ as u;„(A/') — -|-oo for all /3 > 0, as it 
was the case in the previous sections. However, our method does not provide this 
result. 

Proof. Let w be a quasi- free state with Tr(7) < +00. According to [4| Eq. (2b. 26)], 
there exists a Bogoliubov map V:^®.ft— i.e. a unitary operator of the 
form 

■ U J*VJ* 
V JUJ* 



(2.10) V = 



with Tt(V*V) < +00, such that 

(2.11) VT(j,a)V* =T{D,0) 
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with D = diag(Ai)i^o, < +00, and ^ Ai ^ 1/2 for all i. Let uj' be the 

unique HF state associated with r(£', 0) given by Proposition 12.11 Let us also 
denote by U the unitary operator on T lifting V [5^, Theorem 9.5]. Then 

Lui\n){n\) = oj'mn){n\v*) = oj'{\nv){nv\), 

with |rJu) = U|r2). We now estimate |r2u)(f^u| by e^^-^ for any /3 ^ 0. Thus, let 
us fix /? > 0. At the end, we will optimize over f3. By [29, Eq. (67)] we can write 

-1 



i=-K 



i>0 



where (?7i)iez is an orthonormal basis in Sj, {ri.i)-K!^isi-i are eigenvectors of V*V 
for the eigenvalue 1, af + /3f = 1 for all i ^ 0, and {/3f)i^o are the eigenvalues of 
V*V strictly between and 1, which are all of multiplicity 2. We interpret this 
equality by saying that 

where = C?/;, ItJj^) = t], for sC i ^ -1, Sji = C?72i+i ® €7721, iV^i) = ai-A??2iA 
772i+i for i > 0, and = Cry^, jf/'i) = f^i (the vacuum in T{Sji)) ii i < —K. Recall 
that the operator T is the unitary transformation between T{(Bi^i) and 0iJ-'{Sji) 
defined in Appendix \X\ . If —K ^ i ^ —1, then F{S)i) = C ® Cr?i and the matrix of 
the basis (1,77^) can be dominated by 




1 



1 

e-P 



where as usual Afi is the number operator on J^(ioi). If i < —K then F{S)i) = 
C © C?7i and we have 



1 




1 

e-P 



Now if i ^ then J-{S^i) = C Cry2i © C?72i+i © ^7724 A 7721+1- In the basis 
(1, ?72i, V2i+i,mi A 772j+i) the matrix of \ipi){'ipi\ can be dominated by 



IV'z)(V'* 



since with 



/ af 





V 



A 






/3? 



PI 



( 1 








' \ 
























V 











B 



1 





we have B-^'"^ AB^^/"^ (1 + (e^/^ - l)/3f)Idc2 e(^"'-i)^hdc2 . We know that 



MV*V)= i + E2/5f. 
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SO that 



^ T 




T 



i>0 



■ Tr(y y)g-;3AA 



Here we have used that for any operators A,B,C,D, such that ^ A ^ B and 
^ C ^ r> it holds A(S)C ^ B (S) D. We obtain the estimate 



(2.12) sC e 



l)Tr(r>) 



where we have used that uj'{e~^^) ^ e^"^ (M) ^ which is a consequence of the 

equaUty (I2.4p apphed to the HF state uj' . Unfortunately, the estimate (j2.12p is 
not good enough because the constant {e^^ — l)/2 in front of Tt{V*V) is positive, 
while Ti{V*V) represents the number of particles of |riij). We will thus get another 
estimate by exchanging the roles of uj' and The idea is to see |17u) as the state 

and Ld' as the observable. Since ker(D— 1) = {0}, it is well-known that u' is a normal 
state with density matrix G = 2'-iT(M) where M = Z = Tr^(T(A/)), and 

(2.13) T(M) := ®N^oM'^^ : T ^ T 

(see for instance [19l Proposition 6.6 (1)]). Hence we can write 

u:{\Vl){^\) = J{\<^^){^^\) = Tr(Z-iT(M)|l]u)(f^u|) - c^u(^"'T(Af)), 

where wu is the pure state associated with the vector |f2u)- We know that a;u(A/') — 
Tr(y*F). Hence if we can dominate Z-'^T{M) by e"^'-^ for a certain /?' ^ 0, we 
will get another estimate on CLi(|i7) (r2|) by applying equality (|2.6p to the pure state 
wu- Let (/ii)i the eigenvalues of Af, [li = , \ . The spectrum of T(Af) on S^^ is 



a(T(Af)i^,.)=|n^- 



/ C N, #/ = fc 



SO that for / C N of cardinal k we have 

-=[](i-Aon^- 

iel i^O ' iel i^O iel ' i^I iei 

Using ^ Ai 1/2 for all i, we finally get Z-^T{M)iyjk ^ 1/2'= so that 

--2l-zlu;„{Ar) foj. all ^ 0. Choosing 



Equality (12. 6p now implies that aJu(e ■'^) ^ e" 
/?' = In 2, one gets 



(2.14) 



w(|17)(17|) <a;u(e-('"2)^) 



-|Tr(y*y) 
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Interpolating the inequalities (|2.12|) and (|2.14|) we get 

for all /3 ^ and 0^6'^!. We choose such that the coefficients before Tr(Z?) 
and Tr(T^*y) are equal since we have Ty{V*V) + Tr(D) ^ Ty{U*DU) + Tr(T/*(l - 
D)V) = Tr(7) by (piTj) . using UU* sC 1 and 1 - Z? 1. We thus choose 



7 - %e-P + 462/3 ' 

and we obtain 

for all /3 ^ 0. Optimizing the coefficient before Tr(7), we get the desired estimate 
with /3 ~ 0.36443 and ~ 0.308194. □ 

2.6. Bogoliubov-Dirac-Fock case. In this section wc introduce the correct setup 
for studying electron-positron pair creation. Let be a Hilbert space and 11 be 
an orthogonal projection on S^. We also need an anti- unitary operator J as in 
Section 12. 1[ with the additional assumptions that ^ = (i.e. J maps 9) to io) 
and that JIIJ* = 11 or JIIJ* = 1 — 11. The particle/anti-particle spaces are given 
by S)+ = {1- n)^ and = JUS). Notice that jo_ = USj or S)^ ^ [1 ~ X\)S). 
In the context given by Section FOl we have S) = i3A, II = P- and J = i/3a2'rf 
the charge conjugation operator on j^A (i.e. such that J(Z?" + V)J* = —(13° — V) 
for any scalar potential V), with ^ the complex conjugation on j^A- With this 
choice, vectors of Sj^ are interpreted as states with a positive energy relatively to 
the Hamiltonian with an opposite charge. Hence, they represent positronic states. 
Notice that this specific J verifies JP^J* = P° = 1 — P°. In the sequel, we will 
keep a triplet (S), U, J) satisfying the assumptions given above. 

The mathematical description of Bogoliubov-Dirac-Fock states is a special case 
of the well known Araki-Wyss representation [2] (see [HI Section 6.4] for a review). 
Let Jb = J^{^J+) J^i^J-)- For / e Sj+ and g g Sj- we denote by a^(/) and a*_{g) 
the usual creation operators on J^(S)-i-) and J^{S)-), respectively. We now define 
the "creation operator" on Jx) for all / G io by 

rif) ■■= <((! - n)/) (g> Id^(^_) + T (-Id^(^+)) ® a-(Jn/), 

where T is the operation defined by Equation (|2.13p . The operators {ip*{f))f are 
not exactly the usual creation operators in the full Fock space since they create a 
particle in the state (1 — 11)/, and at the same time they annihilate a anti-particle in 
the state JII/, according to the "particle-hole" picture of Dirac's theory. However, 
they still satisfy the CAR thanks to the "twist" T (— IdjF(i5+)) second term. 

A Bogoliubov-Dirac-Fock (BDF) state is, by definition, a quasi-free state uj on the 
C*-algebra 2lo C B{Fo) generated by the {ip{f)) fef) ■ We define the normal ordering 
■i^*{fH{9)- of the operator '4>*{f )^{g) by 

■■rUma)- = <(/+)a+(.9+) ® Id^(^_) -f a;(/+)T (-Id^(^+)) ® al( Jg-) 

+ T (-Id^(y-,^)) a+{g+) (g) a^{Jf-) - ld^(f,+) ® aL{Jg^)a-{J f-), 

where h+ = (1 — H)/i and — lih for all h. It corresponds to moving all the 
creation operators a* to the left of annihilation operators a. For any BDF state w, 
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we define its renormalized one-particle density matrix Q : Sj hy 

{g,Qf)^u;{:rimi9y.)- 
and its pairing matrix p : 9j hy the usual formula 

{pJf,g)=u;{rU)r{9))- 

Recall that we have already defined the 1-pdm 7 by (5, 7/) = Lj{^jj* {f)^{g)). There- 
fore, we have the relation 7 = IT + Q. If A/" is the number operator on To, 

(2.15) TV = ^a\{Lp,^+)a+{ip^^+) ®ldjr(^_) -I- Id^(j--,^) ® a*_{if,^^)a^{if,^^), 

i 

where {Lpi^^)i are orthonormal basis for, respectively, i)^ and i)-, then 

c.(AA)=Tr(Q++-g__), 

where Q++ := (1 — n)(5(l — 11) and Q = IIQII. Hence, uj{N) < +00 is equivalent 

to having Q S &iji{S)) := {Q G 62(-^), Q++,Q S 6i(^)}. Notice also that as 

in the HFB case we have pp* ^ 7 — 7^ = Q++ — Q — Q^, thus p € 62(-^) as soon 

as w(A/') < +00. 

Remark 2.9. The other natural number operator 'Yl,i'^*{fi)'^{fi) gives the total 
number of particles in the system, that is also those of the vacuum IT. If dim(f)_) ~ 
+00, this number is just -|-cx3. However, we only want to count the number of 
particles relative to the vacuum H. That is why the operator M is chosen here: It 
counts the number of "real" electrons and the number of "holes" Lpi- in the 
vacuum. 

The following proposition in an easy adaptation of the arguments given in [3J 
pp. 449-450]. We give the complete proof here to clarify the link between BDF 
states and HFB states. 

Proposition 2.10 (Link between HFB and BDF States). Let {Q,p) G &i^n{^) x 
&2{^) such that 

o^r(H + g,p) = (n + Q ^ . ^(n'+Q)jO ^ ' 

as an operator on S) ® Then there exists a unique, normal, BDF state on 
having Q as renormalized 1-pdm and p as pairing matrix. 

Proof. We are going to construct an HFB state on 2t[io+ ® Sj^] using Proposition 
12. 1[ and then transform it into a BDF state having the desired property via the 
unitary transformation 

T : J^{^+ ® .5-) -> .Fo = J"(io+) (E) J'i^-), 

defined linearly by its action on each iV-particle space 

T(Pt,,+ A • • • A ^,;,,+ A ipj,^- A • • • A = 

A • • • A (p,,, + ) (g) A • • • A <y5j„_,,,-) , 

with the convention Tfl — fi+ (g) Since T leaves the scalar product invariant 
and maps an orthonormal basis for J-{Sj+ ® io-) onto an orthonormal basis for J^q, 
it naturally extends to a unique unitary operator on J-{^+ ® iO-)- It induces the 
following transformation of the CAR: 

(2.16) Ta{ip+®ip_)T*=ilj{ip+®ip_). 
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Define p++ = (1 — n)p(l — 11), p-{ = (1 — n)pn, etc. Suppose also that we are 

in the case where JUJ* = 1-11, and write i5+ = (1 — n)f) = ^3+, i^'L = liSj- We 
now introduce 



7o 



Q++ P++ 
P*++ -JQ — J* 



ao ■■= 



P+- Q+ 
-JQ-+J* p*__ 



where 70 : i3+ ® io+ ^ ^+ ® io+ and ao : Sj°_ ® ^)"_ ^ ^3+ © ■5+- Let us define 



J 



J 
J* " 



and the operator To on {Sj^ ® ^°^) ® by 



To = r(7o, ao) = 



70 



ao 

J*loJ 



We now show that there exists a HFB state utQ on 2t[jo']_ ©i^^].] having To as density 
matrix. Hence we prove that ^ To ^ 1. Let us first write the block decomposition 
of r(n + Q, p) as an operator on Sj^ ® ^)°_ ® ^% ® = Sj ® 



( 



r(n + g,p) 



Q-+ 
p*++ 



Q+- 
i + Q- 

P-+ 
p— 



P++ 

P-+ 
-JQ—J* 
-JQ+-J* 



P+- 
P — 
-JQ-+J* 
1 - JQ++J* 



where we have used that JSj^. — to write the lower right block. Let us consider 
the unitary operator : .Q'j. ©.Q'j. ©.Q° ©io*^ -J> ^^J_ ©^° ©.Q'}. ©.^° whose matrix 
is 



W = 



/I 





V 




1 

/ 



Then, as it was noticed in [3] , we have the relation 

ro^w*r{u + Q,p)w. 

By assumption, we have ^ r(n + Q,p) ^ 1, so that ^ Fq ^ 1 as well. Now we 

have Tr(7o) = Tr((3++ - Q ) < 00 hence by PropositionOwith = io^ ©i3+ = 

© io_ and A = © S^^, there exists a unique, normal, HFB state ujq on 
©io_) with finite number of particles having Fo as density matrix. We define 
a state on Tq via the unitary operator T by 

uj{A) = ujn{T*AT), \fA e BiJ'o)- 



By p.lSp . uj is the quasi-free state with renormalized 1-pdm Q and pairing matrix 
p. Furthermore, uj is obviously normal since luq is normal. This concludes the proof 
of Proposition 12. lOi in the case where JUJ* = 1 — H. If JUJ* = H, the proof is 



the same with VF : i^^ © ^3. 



/I 




V 



^5+ © ij" © i3+ © Sj^_ whose matrix is 



W = 



1 


/ 
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which is the same as in ^ where J was the complex conjugation. With this choice 
of W, W*T{Il + Q,p)W is the density matrix of a HFB state with = -$5+ ® = 
J?. □ 

Corollary 2.11 (Vacuum Expectation of a BDF State). Let uj be a BDF state 

with renormalized one-particle density matrix Q such that Tr((3_|__(_ — Q ) < +oo. 

Then we have the estimate 



(2.17) 



w(|f^o)(f^ol)^e-'^^'-(Q++~ 



where a is the same constant as in Proposition \2.7\ and flo is the vacuum in J-q . 

Proof. Use Proposition 12 . 71 to the HFB mixed state ojq constructed in the proof of 
Proposition [231 □ 

Remark 2.12. In \T7\, the rBDF energy (|1.3p is derived by evaluating the QED 
Hamiltonian on BDF states. As usual in mean-field theories, this energy only de- 
pends on (Q,p). However, since the interaction between the particules is repulsive, 
any minimizer of the rBDF energy has p = 0. This explains why we can take p = 
for the polarized vacuum in the proof of Theorem [TJ 

3. ASYMPTOTICS OF THE POLARIZED VACUUM ENERGY IN STRONG EXTERNAL 

FIELDS 

3.1. Main Results. In this section, we study the asymptotics of the reduced BDF 
ground state energy E{Z) for large Z, where we recall that E{Z) — inf{f j^Qp(Q), Q e 
/C}, with JC defined by Equation (|1.4p . 

Proposition 3.1 (Upper Bound). Let a > and A > 0. Let ly e C n L^{R^,R) be 
such that J l^l — 1 and q — J v 0. Then 



(3.1) 



limsup^l^ -CiaA\qf^^, 



where ci 2-23/331/3,^-4/3 ^ 0.001543 . . .. 

To estimate the convergence rate of E{Z)Z^^/^ ^ we need a further assumption 
on the decay of v at infinity. 

Proposition 3.2 (Convergence Rate). Let a > and A > 0. Letv € CnL^(K^,K) 
be such that J \x\\i>{x) \ dx < +oo, J = 1 and q = J =/= 0. Then there exists a 
constant Z\ = Zi(A,a, \q\, |j|a;|i/|li) > such that 

(3.2) yZ > Zi, E{Z) < -yaA|9|5/3z5/3, 

where ci is defined in Provosition \3. 11 

Remark 3.3. The constant Zi behaves as A3 when A — > oo, as a^'^/^ when a ^ 
and as \q\^^ when g — > 0. It is probably not optimal. 

We also give a lower bound for E{Z), proving that the power Z^/^ is optimal. 

Proposition 3.4 (Lower Bound). Let a > Q, K > 0, and ly e C n L''-{R^,R) with 
J \v\ = 1. Then for all Z > we have 



(3.3) 



E{Z) > -C2aAZ 



5/3 
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where C2 := 2-4/332/3^-1/3 ^ 0.563626 . . .. 

From the asymptotics of E(Z) in Proposition 13.21 we can now derive a lower 
bound on the total number of particles and anti-particles in the polarized vacuum. 

Corollary 3.5 (Average Particle Number of the Polarized Vacuum). Let a > 0, 

A > 0. Leti^ eCnL^(R^,R) with J = 1, / |a;||i^(a;)| dx < +00, and q = Ji^^O. 
Then for any minimizer Q for E{Z) and for all Z > Zi, we have 

(3.4) Tr(Q++ - Q__) ^ CZ''I\ 

I and C is a constant independent of Z , 



where Zi is the same as in Proposition 
given in Equation p.9p below. 

3.2. Proof of the lower bound. We first give the proof of the lower bound in 
Proposition [3]4l which is easier than the upper bound. 

Lemma 3.6. For any Q e K, pq e L°°{M^) and ||/Oq||l~ ^ 

Proof. The proof will only use the fact that any Q G /C is a bounded operator on 
^A- Hence, let Q e /C and V G ^^(R^) n L'^(S?'). Since pq € L^^M^) ^^^^^ ^Yiai 



PQV^Tv^iQV), 



where in the trace V is seen as a multiplication operator on L^(M3,C''). Let us 
denote by IIa the multiplication operator by 1b(o,a) in Fourier space. Since Q is 
an operator on i^A, we have HaQIIa = Q. Now assume that V ^ 0. Then 



PqV 



= |Tro(gnAFnA)| 

^ WUAVU^We, = llVynAlll, ^ J V X ^ttA' 



where we used ||Q|| ^ 1 and the Kato-Seiler-Simon inequality: 

Vp > 2, ||/(x)g(p)||e, ^ {2n)-'^P\\f\\L.\\9\\L.. 

Now if V is not necessarily non-negative, we split V ^ V-^- — V- with V+ — max(y, 0) 
and V- — max(— y, 0). Then we apply the previous bound twice to obtain 



pqV 



PQy+ 



PqV- 



A3 
67r2 



By the density of D L"^ in and by the fact that [L^)' ~ L°°, we get the 
result. □ 

Lemma|32]is crucial to understand the Z^/'^ behaviour of E{Z). Indeed, an easy 
lower bound to fj^op ^rf5DF('3) ^ ~aZ'^D{v, v)/2 for all Q, using the positivity 
of the kinetic energy and completing the square in the other terms. One may think 
that this lower bound would be attained by a Q such that pq ~ Zv, i.e. by a 
state which density of charge compensates the external field. However, Lemma [3^ 
implies that such a state cannot exist in /C, precisely because of the cut-off A. In 
other words, the vacuum cannot "follow" the external field when the field is too 
strong. 
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Proof of Proposition \3.4\ ■ Let Q G /C, then 

ffBDF(Q) > S^'iPQ) ~aDipQ,Z,.) + ^D{pq,pq). 
We also know that for all Q £ K., one has HpqUl^o < ^ —'■ S, so that 

E{Z) ^E{Z,S) :=inf p e C D {R^),\\p\\l^ s$ (J} . 

The variational problem E{Z, S) has the scaling property 

E{Z,6) = Z^E{1,6/Z). 
Define e = 5/Z. We will now show that 

^(l,£)^-J(87r£)V3, 

which then implies Proposition l3.4l Let p be a trial state for E{1, e). Let R > and 
write I • |-i = + 14o with Vi • G Li(K') and • |-i G 

On the one hand 



p{iy-kVi) 



where in the last inequality we used Young's inequality and / — 1. On the other 
hand, 

P(- K.) ^ / Pi.) ( I '-^^y] d. = / Md,, 

where p -.— v-k p and — v{—x) for all x G M"^. Since J|y|>^ = ^(^7 /) with 
/ = (47r)^^5|3,|^fl., as in 21, Proof of Theorem IL3], we have 

P^dy + \d{~p, p) ^ ~\D{f, /) - 

\y\^R \v\ 2 2 2i? 

Notice that 

Drp,p)^M2^rf ^ffiad.^i.(p,p), 

since |P(fc)| ^ (27r)-3/2||,^||^, = (27r)-3/2 for all fc. Hence, 

E{l,e) ^ -2TTaeR^-^. 

ZH 

Optimizing over R, one gets the result. This finishes the proof of Proposition 
lOl □ 

3.3. Proof of the upper bound. Both Propositions 13.11 and 13.21 are proved via 
appropriate trial states. Before turning to the estimates, we start by explaining 
our choice of the trial states. In the sequel, we will assume that J v > 0. The case 
/ < is treated in the same fashion, except for the choice of the trial state (see 
Remark l3.8p . We define the operator Q on by its kernel in Fourier space. For 
ah p,qeR^, let 

Qip,q)^^{p,q)X{p)X{qr, 

with 

7b, q) = (2^)-3/2 f aigAp - 1)F;(p - q)gr{q - 

JB(0,A/2) 
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and 
with 



X{p) = U{p) (1 0)* 



U (p) = a+ (p) - a- ip)p^^, a± (p) = 



'1± 



1 + bl 



2 ■ 



The operator U (p) is unitary on for all p and it diagonalizes D^{p) := a.-p + (3 
as U{p)D'^{p)U{p)* = i/r + ]pF/3- In the definition of 7, we choose < r < A/2 
(a small number which will eventually tend to zero), g,. = r~^^'^g{- /r) with g € 
L^{R^,R), Jg^ = 1 and supp(5) C 5(0,1). We also choose = F{r-) with 
F e i^(K3,M) such that < F{x) < 1 for all x. 

Lemma 3.7. The operator Q belongs to the variational set K. 

Proof. First, notice that Q defines a Hilbert-Schmidt operator on i^(]R^,C^) since 
for instance 



IIQ(T)llW3,C4)<(27r)-=^||5.||i 

^ (27r)-^||5r||i. / 

Jb{o, 



B(0,A/2; 



Al [ dp(\Fr\^i.\ngrf) (p) 



d£ 



A/2) 



\Fr\^^\ngrf 



^ vol(i3(0,A/2)) „ 4 2 



by Young's inequality. It is self-adjoint because Q{p,q) = Q{q,p) for all p,q, and 
we have suppQ(-,-) C (i?(0,A/2) -l-supp^j.)^ C B{0,A)^ hence Q is an operator 
with range in Sj/^. Since for all (fi G Sj/^ we have 



R3 



l{p,<l) {Xiq),(p{q))c4 dq X{p), 



we conclude Q = Qj^ = Q-+ = so that Q = Let ip e S)\. Then 



(<5v, V')l2(m3,c*) = (27r) 



-3/2 



d£((T,5)X^,F*((r,5)X^) 

B(0,A/2) \ / I,2(M3,C*) 



Jb(o,a/2) Jr^ 



dxF{x)\\j'-'iiTeg)X0){x)\ 



S(0,A/2) 



d£ / dp||ff(p-f)X(p)^(p)||^. 



^ II'P||l2(R3,C4); 

where we have denoted by J^~^ the inverse Fourier transform. Hence — P° ^ Q < 

1 - PO. Finally, Tr(0) = f^,jip,p)\\X{p)\\l,dp ^ {2w)-^ J F < +00, so that 
QG6i(i3A)c6i,po(i3A). □ 

Since for any R & K. wc have the formula 
1 



Trc4(i?(p + k,p)) dp, Vfc e 1R^ 



(27r)3/2 

the density of Q can be written as pq = pi+ p2 where for all k, 
pi{k) = {2n)-^VAZ{k)gr*9r{-k) 
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and 

pl(fc) - (2^)-3yA^;(fc) / 9r[p)9r{p + k) {X{p), (X{p + k) - X{p)))^, dp 
with Va := vol(B(0, A/2)) and g := g{--). 

Proof of Proposition \3.1l We start by estimating the terms giving the 
haviour. We have 



(27r)3 7b(o,2A) ' ' 
and 



(2^)V ^(0,2A)' ' " 

We choose r such that (2^)-3/2gZ ^i^r-^ = (^)%-5, ie = (^) 
The constraint r < A/2 is equivalent to qZ > (27r)3/^ IT' which is automatically 
satisfied in the limit Z — > +00. We will come back to it in the proof of Proposition 
We thus get by the dominated convergence theorem 



(3.5) lim (^aD{pi,Zi^) + -D{pi,pi)) = 2-''/^3-'/\''^^^aAq'^^x 

Z— J-+00 \ 2 / 



F(% * g{-k) ^j^^l f |f(%*g(-fc)P 



■(0,2A) I^P 2 7s(0,2A) 1^1 / 

We now want to optimize the right side with respect to g and F. We choose 
F = lB(o,a) and g = (3(47r)"i6"3)i/2lB(o,6) with a > and 6 € (0,1), and we 
optimize over a and 5. Since g*g{k) — {vol{B{0,b)))^^vol{B{0,b) n B{~k,b)) for 
all fc and since B{-k/2, b - |fc|/2) C 5(0, b) n S(-fc, 6) for aU |fc| < 26, we have 



g*gik) ^ lB(o,2b)ik) [I 25 



For all |fc| ^ 26, we also have 



Therefore, 



s(o,2A) I'*''!" v27r V3 / JB(ti,2b) \ 26 

V27ra^6 ( - - ab 



3 

Then, using |F(A;)| ^ 47ra3/(3(27r)3/2) and \gi<g{k)\ < 1 for all fc, we obtain 



(3.6) / 
Je 



|F(%^g(-fc)P ^ 16^6^ 

i3(0,2A) ^ 9 
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Hence for the Z^l'^ term we find that 

Z— i-+oo \ 2 / 



2\ 8 



af'h V27r \ ah - - , . 
V V 37 9 

Opthnizmg the right side over aU a > and 6 £ (0, 1) we get the result because 
min f ^2^ f a6 - - V ^ _2-'5/6^'/3^5/6 ^ 

a>0 V V 3 / 9 / 

6e(o,i) \ \ / / 

It now remains to prove that the other terms in £'j^'^p((5) are of lower order 
than Z'"^/-^. We begin with the kinetic energy, 



Trop°Q) v/r+A^Tr(Q) = {27r)-\-^ I F = {27r)-' a' qZ. 



A3 

For the terms involving p2, we first use that for all p, k, 

\\Xip+k)-X{p)\\ ^ \\U{p+k)-U{p)\\^^\k\. 



Consequently, for all fc. 



|P^2(fc)| -^^r-'\F{k/r)\\g*g{~k/r)\\k\. 
Using this bound together with the estimates leading to p.6p . one finds that 

\D{p2,Zl,)\ ^ /3iA2gl/3^4/3^ 
\D{puP2)\ ^ iS^A^q^/^Z^/^ 

\D{p2,P2)\ ^ Pz^-\Z, 

/3i = 2i9/33-2/377r23/6a3ft2^ 

/32 = 2"/63-8/37^-7/3^652^ 

/?3 = 2-3/23-4727r-5/2a6^3_ 
limsup^^ lim ^i%M = -ciaA^s/s. 

Z^+^ Z5/3 ^ Z^+oo Z5/3 



Remark 3.8. In the case ^ v < Q, the proof is the same except that we take 



with 



We conclude that 



(0 1) 



□ 



X{p) - U{p) 
and 

QiP. q) = -lip, q)X{p)X{q)\ yp, q G R\ 
The trial state Q now verifies Q — Q , so that pq is locally negative around 0. 

To estimate the convergence rate of E{Z) towards — ciaA|g|^/3^5/3^ ^^y[ use 
the first moment / IxWiy] to control the convergence of to t?(0). 
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Proof of Proposition \3.2l We assume q > (the case q < follows from obvious 
modifications). We split the term D{pi^Zv) into 



(2^)3/2 7b(o,2A) \k\ 

with 



JB(0,2A) \'^\ 

We use that \v{k) - 9(0)| < (27r)-3/2|||2,|j^||^|;j| fo^ ^ estimate 

«C 25/33-5/3^-ll/6„3^2|j|^|^|j^y^2|^|l/3^4/3_ 



47r 



Hence we have for all \q\Z > ^^^^3^^ -g-, 

(3.7) i?(Z) ffBDF(Q) < -ciaA|g|5/3z5/3(l - BiZ-i/'^ - B2Z-^/'), 

with 

Si = 228/33-27r-i/V62|||x|i.||iA|9|-4/3 + 24/33-17^31/6^3^2^1^1-4/3 

+ 251/63-37^-1^6^2^1^1-1/3 

and 



B2 = 237/63-1/3^4/3^3 VI +^ |^|_2/3 ^ 231/63-ll/372^-7/6a6;,3^2|^|-2/3^ 

We furthermore have 1 - BiX - B2X'^ ^ 1/2 for X s$ ATq with 



therefore Proposition 13.21 holds with 
(3.8) Zi := max (x-^, ^ ^ 



(2^)3/2 3|g|j 



□ 



Proof of Corollarv \3.5\ For any Q e /C, we have the estimate [HI Lemma 1] 

^(pq,Pq)'/'^C'a|1Q||i,po, 

with Ca < +00. From the constraint — ^ Q ^ 1 — P° we also get Q++ — Q ^ 

Q^, hence defining X :— [Tr(Q++ — Q )]^^2 we have 

2II0+-III, = 2II0-+III, - IIO+-III, + IIQIII, x^ 

hence 1 1 Q 1 1 1 po ^ X^ + ^/2X . Using the Cauchy-Schwarz inequality for the Coulomb 
scalar product we also get that for all Q G /C 

D{pQ,Ziy) ^ CaZMc{X^ + V2X). 

From this estimate and the inequality p.2p . we can see that for all Z > Z\ and for 
all minimizer Q™'^ for E{Z), 

-aC^ZMciX^ + V2X) ^ £fBDF(Q™') = E{Z) ^ -yaA|g|5/3z5/3. 
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Hence, 

X ^ ^ (Vl + 2a- 1) 
v2 

witha- §i^£^Z2/3 so that 



where 




using that x (vT~Pa; — 1)/Va; is increasing on (0, +oo). □ 

4. Proof of Theorem [T] 

Proof of Theorem]^ Let Q be any minimizer for E{Z), and let uj be the unique 
BDF state on J^{S)+) ® having Q as its generahzed 1-pdm and p = as its 

pairing matrix, defined by Proposition 12 . f 01 Then by Corollarv l2.f f 1 we have 

Now by Corollary [231 we know that for all Z > Zi, 

Tr(Q++-Q— ) 

Thus pz^l- e"'"^''' with k := aC. □ 

Appendix A. Product States 

Given a (at most) countable family of separable Hilbert spaces (jOi)^^^ and a 
family (uji)i^fq of quasi- free states such that is a state on J-{S^i) for all i e N, 
we want to give a meaning to the product state ^i^noJi as a quasi-free state on 
•^(ffiiGN-^i) • We first consider the unitary transformation 

AjeJ., VjM A • • • A Aje,/.^ 'P].^l. ^ A^g j.^ VjM «) • • • ® Aje j., Vj',**. . 

where ii < ■ ■ ■ < ik are elements of N and for each 1 ^ ^ A:, is 
an orthonormal basis for Jo^^ and Ji^ C N is finite. We recall the definition of a 
product state on a tensor product of C*-algebras [Til Proposition 2.9]. 

Theorem 2. Let (2li)igN be a collection of (unital) C* -algebras and let uji be a 
state on 21^ for all i. There exists a unique state on ^iSti, denoted by ®iU)i such 
that for any {Ai-^ , . . • , Ai^) G 21^^ x • • • x 21^^ with any indices ii < ■ ■ ■ < ii, we have 
(S)iUJi {Ai^ (K) . . . (g) J = Ilk=i ^tki^ik)- 

Remark A.l. Recall that the tensor product 0^21^ is defined as the inductive limit 
of the C*-algebras j2ti with J finite, or equivalently as the completion (for a 
certain *-norm) of seen as a tensor product of unital algebras 

In the particular case where 21^ — 2l[.$3i] for all « e N, the tensor product 0^21^ is 
also a CAR algebra generated by the operators 

S(/,) := T (-Id^J ® • • • ® T (-Id^^._ J ® a,(/,), 
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where j G N, fj E Sjj and aj is the annihilation operator on J-{Sjj) for aU j. Notice 
the "twisting" operator T (—Id) on the left used to ensure that the (a(/j)) f^j satisfy 
the CAR. Notice also that we have the relation 

where a{fj) is the usual annihilation operator on T{(Bi^i)- 

Proposition A. 2. Let (cji)^^^ be a collection of quasi-free states such that Ui is 
defined on ^[Sji] for alii. Let{'^i,ai)i he the collection of their density matrices, and 
assume that X]iTrjji(7i) < +00. Let u he the unique quasi-free state on 2l[®ii5i] 
having 7 :— Oiji as 1-pdm and a := OiOi as pairing matrix. Then we have 

VAe 2t[e,^),], oj{A) = (E),oj^{TAT*). 

As a consequence, the product state ®iLOi can he considered as a state on 21[®ii3i]. 

Remark A. 3. The state lo is well-defined by Proposition l2.11 since Tr(7) = Tr^^ (7^) 
is finite. 

Proof. We will show that the state ^iOJi is quasi-free on the CAR generated by 
the {a{fj))f-j. This proves the result because the state ®iOJi{T ■ T*) is then quasi- 
free and one easily shows that it has ©i7i as 1-pdm and ©iCti as pairing matrix. 
Since the density matrices determine uniquely the quasi-free state, we must have 
®iijJi{T ■ T*) = uj. Therefore we only have to show the Wick relation for the state 
®iUii. Let us first notice that it is enough to prove it for products of the form 

a« (/i-))...a« (/ir;), 

with ii<--- <iN,ktEn, f^'"^ e iOj, for aU 1 s$ p 2ki, 1 ^ s$ iV, and where ji 
means star or no star. This means two things: 

(1) We can restrict to ordered products with respect to the decomposition 
fflieN-^j: The 2fci first creation/annihilation operators a" (^fp^^^ with 1 ^ 
p ^ 2fci all create/ annihilate particles belonging to iOi^, then the follow- 
ing 2fc2 creation/annihilation a" (^fp^^^ with 1 ^ p ^ 2/c2 create/annihilate 

particles belonging to S^i^ , etc. We can always order in this way any product 
of a's because Wick's formula does not depend on the choice of an ordering. 

(2) We can restrict to products where there is an even number of particles 
created/annihilated in each space Sji^,. If not the case, both sides of Wick 
relation can easily be shown to vanish. 
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This being said, let us compute the left-hand side of Wick's formula: 

X = [a« . . . a« . . . a« (/f • • • a» {ftj)' 

= ®iu:, [t (-Id^„)^"-2'=^ ® • • • ® T (-Id^,.^_,)^"=^''^ 
®a« (/^') • • • a« (/i-') T (-Id^ J^"- ''^^ • • • 



= n E (-1) 



e(7rf ) 



Let us reindex the product 



(/^'/(l)) "'^ (-^^^(2)) 

(/^/(2/c<-l)) 



X w,: 



(if) 



S« (/i) • • - S* (/2fcJ a« (/2fe,+i) • • • a« (/2(fc,+...+fe,)) , 
such that the sum of Wick's formula (with K := ki + . . . + fcAr), 
y = ^ [a« (/,(!)) a« (/,(2))] • • • {^iUJi) [a« (/.(2if-i)) S« (/.(2k))] 

contains only the terms with a permutation tt leaving invariant every interval of 
the form [2ki + 1, 2A:^+i], corresponding to the Hilbert space Sjig, 1 ^ £ ^ N. This 
is proved by induction on K: We introduce 

W^{ei . . . e2K) ■■= ^ (-l)^(''V(e7r(l)e7r(2)) • • • P{en{2K-l)en{2K))- 

For TT € S2K1 we have 7r(l) = 1, hence we have the induction formula 

2K 

(ei . . . e2K) = ^(-l)V(eiei)M//-i(e2 . . . . . . e2K). 



i=2 



In our case where Cj = a'(/j) and p = ®iLOi^ we see that p{e\ei) 7^ if and only if 
i ^ 7k\. Hence, by induction on the only permutations tt S iS2i4: giving rise to 
a non-zero term in Y are those which leave invariant the intervals [2A;^ + l,2fc^+i] 
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with 1 ^ ^ ^ A^. We can thus write 

N 



y = (-1) 



AT 



n E (-1) 



H) a" (/^/(2fe<-l)) (/^/(2/s^)) 

(■^^/(Ifc^-l)) (/l/(2fc«)) ^ ■ ■ ■ 



TT^eiiT, 

• • • X 

= X. 

This proves that the state ^jW, is quasi-free, which concludes the proof. 
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